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Abstract. A dynamically preferred quasi-local definition of gravitational energy is 
given in terms of the Hamiltonian of a '2+2' formulation of general relativity. The energy is 
well-defined for any compact orientable spatial 2-surface, and depends on the fundamental 
forms only. The energy is zero for any surface in flat spacetime, and reduces to the 
Hawking mass in the absence of shear and twist. For asymptotically flat spacetimes, the 
energy tends to the Bondi mass at null infinity and the adm mass at spatial infinity, taking 
the limit along a foliation parametrised by area radius. The energy is calculated for the 
Schwarzschild, Reissner-Nordstrom and Robertson- Walker solutions, and for plane waves 
and colliding plane waves. Energy inequalities are discussed, and for static black holes 
the irreducible mass is obtained on the horizon. Criteria for an adequate definition of 
quasi-local energy are discussed. 

1. Introduction 

It comes as a surprise to many that there is no agreed definition of gravitational energy 
(or mass) in general relativity. In Newtonian gravity, there is a material density which 
may be integrated over a spatial 3-surface to give the mass, and the Poisson equation and 
Gauss theorem then yield an expression for the mass inside a 2-surface, measured on the 
2-surface — a so-called quasi-local mass. Finding a similar quasi-local expression for the 
Newtonian gravitational energy is impossible, since the energy density, unlike the material 
density, is not a total divergence. In general relativity, material mass is merely one aspect of 
the energy-momentum-stress tensor, and gravitational energy, if well-defined at all, is non- 
local, as follows from the equivalence principle. The gravitational field can be measured 
by the geodesic deviation of two observers, but a single observer cannot distinguish it from 
kinematical effects. Equivalently, curvature cannot be measured on a point or line, but 
requires a 2-surface at least. 

In the case of an asymptotically flat spacetime, the Bondi and adm masses are well- 
defined asymptotic quantities which are generally accepted as the total mass of the space- 
time measured on spheres at null infinity and spatial infinity respectively (Bondi et al 1962, 
Sachs 1962, Arnowitt et al 1962). The existence of these asymptotic definitions, and in 
particular the Bondi mass-loss result, indicate that gravitational energy is not meaningless, 
and that there is a real physical process which is responsible for radiating gravitational 
energy to infinity in such a spacetime. Since an appropriate definition of this energy cannot 
be found locally, a quasi-local definition is sought, where 'quasi-local' is here taken in the 
conventional sense of referring to a compact orientable spatial 2-surface, usually spherical. 
This is the natural finite analogue of the asymptotic definitions. 
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Hawking (1968) denned a quasi-local mass which has various desirable properties: it 
is zero for a metric sphere in flat spacetime, it gives the correct mass for the Schwarzschild 
solution on a metric sphere, it tends to the Bondi mass asymptotically in a static asymp- 
totically flat spacetime, and it tends to zero as a sphere in any spacetime is shrunk to a 
point (Eardley 1979, Horowitz and Schmidt 1982). Unfortunately, it is non-zero for generic 
2-surfaces in flat spacetime. This drawback is easily corrected by adding a certain term, 
as is shown subsequently. 

Penrose (1982) emphasised the conceptual importance of quasi- local mass, and sug- 
gested a twistorial construction. This can be properly defined if the 2-surface can be 
transplanted into a conformally flat spacetime (Tod 1983, 1986), but the transplant can- 
not be done generically without damaging the structure of the 2-surface (Heifer 1992). 
Bergqvist (1992a) reviews various other attempts, of which that of Dougan and Mason 
(1991) comes closest to being a well-defined functional of 2-surfaces, being well-defined 
generically, but breaking down for the important case of marginally trapped 2-surfaces. 
Bergqvist also shows that seven different definitions all give different results in two simple 
examples, namely the Reissner- Nordstrom solution and the Kerr horizon. This is an unfor- 
tunate situation which indicates the need for a more critical approach to what constitutes 
an adequate definition. Currently, there seems to be no demand for a quasi-local energy 
to be even well-defined, let alone geometrically or dynamically natural. 

The main purpose of this article is to present, in §3, a dynamically preferred quasi- 
local energy which is essentially the '2+2' Hamiltonian of the Einstein gravitational field. 
This is a well-defined functional of compact orientable spatial 2-surfaces, and is geomet- 
rically natural in the sense of depending only on the fundamental forms of the 2-surface, 
introduced in §2. The energy is shown to vanish for any surface in flat spacetime in §4, 
to tend to the Bondi mass at null infinity in §5, and to tend to the adm mass at spatial 
infinity in §6, if the spacetime is asymptotically flat. Uniqueness on purely geometrical 
grounds is considered in §7, examples in §8, energy inequalities in §9, and criteria for an 
adequate definition in §10. 

2. Geometry of 2-surfaces 

It is appropriate to begin with a review of the geometry of a compact orientable spatial 
2-surface S embedded in spacetime, according to the '2+2' formalism developed by the 
author (Hayward 1993a), which describes null foliations of such surfaces. In this formalism, 
a basis (e" , e^) for S is completed to a spacetime basis (u a , v a , e", e^) by Lie propagation: 
(£ u , C v )(u a , v a , e", e%) = 0, where C denotes the Lie derivative. The commuting vectors 
(u a ,v a ) are referred to as the evolution vectors, since Lie propagation in these directions 
enables the spacetime in a neighbourhood of S to be developed. (This is analogous to 
the Cauchy problem, where the single evolution direction may be decomposed into a lapse 
function and a shift vector.) The choice of evolution vectors is partially fixed by demanding 
that the 3-surfaces (£ U S,£ V S) are null, u a l a = v a n a = 0, where the null normals (l a ,n a ) 
satisfy 

V [a / 6 ] = V [a n b] = 0, l a l a = n a n a = 0, 
l a n a = -e m , h ab l b = h ab n b = 0, g ab = h ab - e~ m (l a n b + n a l b ), 
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and g a (, is the spacetime metric, h a b the induced 2-metric of S and m is called the scaling 
function. The remaining freedom in the evolution vectors is given by the shift 2-vectors 
r a = hfrU b and s a = hiv . The metric then decomposes into 



r c r r c s — e 

(Jab - ! r c s c - e~ m s c s c 



r, 
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in the above basis. The freedom to interchange or boost the null normals is left open. 
The dynamically independent first derivatives of the metric consist of the expansions 

9 — £ u — r h c di 9 2^ ^v—shcdi 

the (traceless) shears 

®ab = h a h^C u — r h C( i 2^ a b^ ^u—r^cd-, o ab = h a h^C v — s h C( i ^h^h C v — S h c d, 
the 'inaffinities' 

v = C u - r m, v = C v - S m, 
and the 'twist', or 'anholonomicity', or commutator of the null normals, 

oo a = \e m h ab (C u s b - C v r b - C r s b ). 

These fields encode the extrinsic curvature of S, or equivalently the momenta conjugate to 
the configuration fields {h a b, r a , s a , m). In Hayward (1993a), the vacuum Einstein equa- 
tions are written in a first-order form in terms of (h a b, r a , s a , m, 6*, a a b, &ab, v, z>, u a ), which 
constitutes a '2+2' analogue of the adm '3+1' formalism. The initial data for the associ- 
ated characteristic initial value problem can be taken as (h a b, m, 9, 9, ui a ) on S, (<J a b, v) on 
C U S and (a a b, z/, s a ) on C V S, with r a being prescribed over the whole 4-d patch. Of this 
data, r a , s a , m, v and v represent coordinate freedom on the respective surfaces, and will 
be set to zero henceforth. 

The remaining dynamical variables are, in traditional language, the first fundamental 
form hai, the (null) second fundamental forms a a b + 9h a b and a a b + 9h a b, and the normal 
fundamental form u a (cf Spivak 1979), with the null normals fixing a preferred basis for 
the fibres of the normal bundle, up to boosts and interchange. These fundamental forms 
are the only geometrical invariants of a 2-surface in spacetime, and are also the dynamically 
independent parts of the gravitational field referred to a 2-surface. It is therefore most 
natural to seek a definition of quasi-local energy in terms of these forms. In fact, of the 
various definitions, this is the case only for the Hawking mass, 



1 / A 
8tt V 16tt J s ' 

where n is the area 2-form and 1Z the Ricci scalar of h a b- 



3 



3. Definition of the energy 

To measure the energy of a field on a spatial 3-surface, the Hamiltonian of a '3+1' 
formulation is usually taken. On a compact orientable spatial 2-surface S, the analogous 
dynamical object for the Einstein field is the Hamiltonian 2-form, which may be quoted 
from Hayward (1993a) as 

8nH = -n(K + 66- \o ah b ah - 2u a u a ), 

where the notation and choice of evolution are as in §2, and the 8tt has been inserted to 
agree with the units G = 1 used in most papers on quasi-local energy. This is a preferred 
dynamical quantity in the sense that variation of the (full) Hamiltonian yields the Einstein 
equations in a form adapted to null foliations of spatial 2-surfaces. The Hamiltonian is 
obtained from a Lagrangian which is the Einstein-Hilbert Lagrangian in '2+2' form, up 
to a total divergence. In a '3+1' context, the energy is defined as the integral of the 
Hamiltonian over the 3-surface, but in a '2+2' context, f s 7~t needs to be multiplied by a 
length to give energy units, and the only natural length scale is given by the area 

A= /V 

Js 

The proposed definition of the gravitational energy of a 2-surface, or its quasi-local energy, 
is then 




where the factor has been chosen for agreement with the Schwarzschild mass. Thus 

E= h{w, j^ + 0~0-\a ab ~a^-2u a u«). 

It should be noted that the expression is valid for any one 2-surface, but that when consid- 
ering the variation of energy between surfaces, the quasi-local coordinate freedom cannot 
be fixed as in §2 and the full Hamiltonian should be used. 

Various immediate observations concerning E can be made. Firstly, the total energy 
of two surfaces S = S\ U S 2 is greater than the sum of the individual energies, since 



A = A 1 + A 2 , 




so that E > Ei + E 2 . The difference can be interpreted as an interaction or binding 
energy, showing that it is energetically favourable for separate concentrations of energy to 
coalesce. Secondly, the term in the Ricci scalar may be integrated by the Gauss-Bonnet 
theorem {e.g. Spivak 1979), 

/ ^ = 8tt(1-<7), 
Js 
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where the topological invariant g is the genus (number of handles) of S: g = for a 
sphere, g = 1 for a torus, etc. Thirdly, E reduces to the Hawking mass Mh in the shear- 
free, twist-free case. Finally, the shear term in E is precisely the addition to the Hawking 
mass required to yield zero in flat spacetime, as follows. 

4. Flat spacetime 

Recall first the contracted Gauss equation (e.g. Spivak 1979), 

1Z + 98 — \o ab d ah = h ac h bd R a bcd, 

which is a purely geometrical equation describing the embedding of S, with R a bcd being 
the Riemann tensor of g ab . 

In flat spacetime, commutativity of the null normals means that the twist vanishes, 
uj a = 0, and the Gauss equation becomes 

K + 09- \o ah ~a ah = 0, 

so that Tt and E vanish. 

Alternatively, note that the expansions and shears can be taken as equal and opposite, 
e = -e, a ab = —a abl by fixing the boost freedom of the normals. For a 2-surface which 
lies entirely in one Euclidean 3-plane, Euclidean surface theory (e.g. Spivak 1979) can 
be applied. Here, the Ricci scalar is twice the Gaussian curvature, and so is defined by 
1Z = 2k + k,-, where the principal curvatures k± are the roots of the eigenvalue equation 

det (y^(v a b + 0h ab ) - Kh ab ^j = 0, 

since the null second fundamental form is \[2 times the Euclidean second fundamental 
form. This yields 

V2k± = 0± \a ab a a \ 

and so 

Tl - 6 2 + \a ah a ah = 0, 

so that H and E vanish. 

Since a ab CF ah < 0, it also follows that the Hawking mass is non-positive in flat space- 
time, Mh < 0, vanishing only for shear-free surfaces. This problem with the Hawking 
mass was well-known, and its solution is simple: add the appropriate term in the shear, 
as determined by the Gauss equation. A similar argument shows that the mass of Geroch 
(1973) is also non-positive in flat spacetime. Horowitz and Schmidt (1982) showed that 
the Geroch mass is always less than the Hawking mass. 

5. Null infinity 

Asymptotically flat spacetimes were first studied in the axisymmetric case by Bondi 
et al (1962), and in general by Sachs (1962). They defined an asymptotic mass Mb on 
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slices of null infinity, and showed that it decreased to the future, which is interpreted as a 
loss of mass due to gravitational radiation. Briefly sketching the approach: on a foliation 
of null surfaces labelled by u, spherical polar coordinates (r, ip) were introduced, such 
that each surface S of constant (u,r) had area form \i = r 2 sin^cfr? A d<p, and hence area 
A = Aixr 2 . Asymptotic expansions in negative powers of r were assumed, and the Bondi 
mass Mb was defined as the integral over sind dfi A dip of a certain coefficient, giving a 
generalisation of the Schwarzschild mass. In the Appendix it is shown that the Bondi mass 
may be written as 

M B = lim ^\[^- [ ^h ac h bd C abcd , 

r^oo g7T V 1D7T J g 

where C a b c d is the Weyl tensor of g a b- Thus the Bondi mass may be described as the average 
asymptotic Coulomb part of the gravitational field. The advantage of such an invariant 
expression in the present context is that the radius r does not appear explicitly except as 
the limit parameter, so that the limit could simply be removed to yield a quasi-local mass 
corresponding to Mb- Using the vacuum Gauss equation 

1Z + 99 — \a a ba ah = h ac h bd C a bcd, 

the Bondi mass may be rewritten in terms of the fundamental forms as 

M B = lim ^\\4r [ v{n + 9~9-\a ah ~a ab ). 

r^oo SlV V 1D7T J s Z 

Comparing with the quasi-local energy E, the asymptotic behaviour of the various terms 
may be translated from the spin-coefficient expressions (e.g. §9.8 of Penrose and Rindler 
1986) as 

n = 2r- 2 + 0(r- 3 ), 99 = -2r~ 2 +0(r~ 3 ), a ab d ab = 0(r" 3 ), u a u a = 0(r" 6 ), 
so that the asymptotic limit of E exists, E = 0(1), and equals the Bondi mass, 

M B = lim E, 

r^oo 

since the twist term tends to zero. Thus E satisfies the criterion of tending to the Bondi 
mass asymptotically, if the latter exists. 

Note that the expression for Bondi-Sachs mass given by Penrose (1963, 1964, 1967) is 
an asymptotic limit of the Hawking mass, 

M B = lim M H , 

which appears to differ from the previous expression by a term in the shears. This is due to 
a different limit being taken, using a foliation based on an affine parameter f = r + 0(r _1 ), 
as is explained in detail in the Appendix. The affine expression for the Bondi mass is 
used in most work on asymptotic flatness (e.g. Newman and Tod 1980, Winicour 1980, 
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Wald 1984, Penrose and Rindler 1986, Christodoulou and Klainerman 1989, Stewart 1990). 
Other expressions are given by Bramson (1975) and Streubel (1978). See also Winicour and 
Tamburino (1965), Tamburino and Winicour (1966), Geroch (1977), Geroch and Winicour 
(1981) and Walker (1982) for the 'linkage' formulation of asymptotic mass and angular 
momentum. 

The positive-mass proofs of Ludvigsen and Vickers (1982), Horowitz and Perry (1982) 
and Reula and Tod (1984) all use the affine version, as can be seen by comparing their 
expressions for mass with the definitions of Ludvigsen and Vickers (1981) and Horowitz 
and Tod (1982). A positive- mass proof has been given by Schoen and Yau (1982) using 
the original approach. 

6. Spatial infinity 

Arnowitt et al (1962) defined an asymptotic mass Madm at spatial infinity, which 
is interpreted as the total mass of an asymptotically flat spacetime. The definition may 
be written in a more coordinate-independent way (Ashtekar and Hansen 1978, Ashtekar 
1980, Ashtekar and Horowitz 1982) as 



M ADM = lim f\/-f / fih ac h bd C abcd , 



p^oo 87T V 16"7T J g ' 

with p denoting the area radius of a family of spheres S which approach spatial infinity in 
a spatial 3-surface (Appendix). As at null infinity, the twist term in E disappears in the 
limit, so that 

Madm = lim E. 

Thus Mb and Madm are null and spatial limits of the same quantity E, or equivalently 
of the average asymptotic Coulomb part of the gravitational field. This agreement has 
apparently not been noticed previously. Indeed, in relating Madm to the limit of Mb, 
Ashtekar and Magnon-Ashtekar (1979) found it necessary to require the vanishing at spatial 
infinity of the shear term that apparently distinguishes the area-radius definition from the 
affine version (Ashtekar and Horowitz 1982). Note that it is not clear whether Madm 
is the limit of the Bondi mass Mb at spatial infinity, since the limits are from different 
directions and may not coincide. The situation is simply that if Mb or Madm exist, then 
they are the appropriate limits of E. 

7. On uniqueness 

The results of §4-6 indicate that the freedom to choose a sensible quasi-local energy, 
even on purely geometrical grounds, is quite limited. More precisely, assume that the 
energy can depend only on the fundamental forms of the 2-surface, including invariance 
under interchange and boosts between the null normals. Assume also that the energy is 
second-degree in g ab , i.e. a linear combination of second derivatives and quadratic first 
derivatives (cf the discussion in §11.2 of Wald 1984). Then the only such functions are 
linear combinations of 1Z, 99, o a b<7 ah and u a u a . The coefficients of 1Z and 99 are determined 
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by consideration of metric spheres in flat spacetime and in the Schwarzschild solution, 
and the coefficient of a a b& ab is determined by demanding that the energy vanish in flat 
spacetime, as in §4. This automatically yields the Bondi mass Mg at null infinity, as in 
§5, and the adm mass at spatial infinity, as in §6. This leaves only the coefficient of uj a u a , 
which is not so clearly determined, but can be justified as a contribution due to angular 
momentum (Hayward 1993c). The above assumptions are hardly compulsory but seem 
quite plausible, and it is sensible at least to exhaust such possibilites before resorting to 
more radical suggestions. It is certainly remarkable that E automatically satisfies these 
desiderata. 



8. Examples 

In calculating the energy for particular spacetimes, it is often most convenient to use 
the full Hamiltonian. Alternatively, in simple cases a transformation may be sought in 
which the line-element, evaluated at the 2-surface S in question, takes the form 

ds 2 = -2d£drj + h ab dx a dx b , 

where (x l ,x 2 ) are coordinates on S and (£,rj) are affine parameters; this is referred to 
subsequently as the standard form. The fundamental forms may then be calculated by 
noting that C u = d/d^ and C v = d/dt] in the basis used. Also useful are the expressions 
[id = /i£ and fi9 = n v . 

In the particular case of spherical symmetry, a sphere of symmetry has zero twist and 
shears, and the same results are found as for the Penrose mass according to Tod (1983) 
or the Hawking mass. Here spherical polar coordinates (r, ip) may be taken such that 
[i = r 2 sin $ dd A d(p, and conseqently A = 4/rr 2 , 1Z = 2r~ 2 , 6 = 2r _1 r^, 9 = 2r~ 1 r v , and 
so 

E = r(r^r v + \). 

This definition of gravitational energy is often regarded as standard in spherical symmetry, 
and was originally given in a different form by Misner and Sharp (1964). A few examples 
follow, with interpretive comments. 

The Schwarzschild black-hole solution 

ds 2 = - (l - dt 2 + (l - ^ 1 dr 2 + r 2 {d$ 2 + sin 2 dip 2 ) 

can be put in standard form by taking 

1 / 2m\ 1/2 1 / 2m\- 1,2 

^ = -^ = -^r 2 { 1 --) ' ^ =t - = T2\ l -~) ' 

and so E = m. It should be emphasised that an interior solution matched to the 
Schwarzschild exterior need not have total material mass m as measured on a spatial 
3-surface. In this sense E is not sensitive to the distribution of matter inside S, being 
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only a single number, but rather measures the effective active gravitational mass as felt 
on S. Another way to see the necessity of considering the effective mass is to note that 
the maximally extended Schwarzschild solution has zero material mass, being a vacuum 
solution. 

The Reissner-Nordstrom charged black-hole solution 

2m e 2 \ ,. 2 / 2m e 2 \ 2 2/r Q 2 , • 2„ a 7 2\ 



ds z = - 1 + -z)dr+[l + -2 dr z + r\d$ z +sin z $d(p z ) 



can be put in standard form by taking 



2\l/2 i / .2N-1/2 



1 / 2m e 2 V^ 1 / 2m e' 



and so E = m — |e 2 r _1 . The same result was obtained for the Penrose mass by Tod 
(1983), who noted that the term in e agrees with the linearised limit. If this is regarded as 
a correction to m which yields the effective mass, then there is the interesting result that 
the field becomes repulsive close to the r = singularity, a result which is also indicated 
by the behaviour of geodesies (§5.5 of Hawking and Ellis 1973). For a non-extreme black 
hole, e 2 < m 2 , such negative E occurs inside the inner horizon, which is widely regarded 
as unstable due to an 'infinite blue-shift' effect, in which case the inner region would not 
exist in practice. Nevertheless, the example indicates that a negative energy could be 
interpreted in a quasi- Newtonian way as a net repulsion of the 2-surface (cf Kulkarni et 
al 1988). 

The Robertson- Walker cosmological solutions are given by 

ds 2 = -dt 2 + a 2 (dif; 2 + f 2 (d$ 2 + sin 2 dp 2 )) , 

where a(t) and 

{sin if; for k = 1, 
V> for k = 0, 

sinhip for k = —1, 

where k labels the spherical, flat and hyperbolic cases respectively. This may be put in 
standard form with radius r = af by taking 

The field equation a 2 = |7ra 2 p— k, where p is the density, then yields E = ^nr 3 p. Recalling 
that A = 4nr 2 , E is just the product of density and the 'area volume' associated with the 
area radius r, rather than the actual volume of a 3-surface inside S. Again, this leads 
to the interpretation of E as the effective mass felt at S. Note that E is independent of 
the pressure, which could for instance be negative, as in the de Sitter and anti de Sitter 
solutions. 
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Although quasi-local mass is usually considered only for spheres, the definition of E 
applies to any compact orientable 2-surface, and it is interesting to compare the results 
obtained for a flat torus in a plane-symmetric spacetime with topological identifications. 
The general plane-symmetric line-element is given by Szekeres (1972) as 

ds 2 = -2e~ K dudv + e~ p (e Q coshWdx 2 - 2 sinh Wdx dy + e~ Q cosh Wdy 2 ), 

where (K, P, Q, W) are functions of (u,v), and the toroidal identifications are 

(x,y) = (x + x ,y + y ). 

The metric may be put in standard form by taking 

«,£ = v v = e K l 2 ', u v = V£ = 0, 

and it is straightforward to calculate fx = e~ p dx Ady, A = e~ p xoyo, 1Z = 0, 99 = e K P u P v , 
a ab a ab = 2e K (Q u Q v cosh 2 W + W U W V ), oo a oo a = 0, and hence 

E = \ (^) 3/2 e K ~ ZPI \PuP v - QuQv cosh 2 W - W U W V ). 

Plane wave spacetimes are given by the above metric depending only on u, or only on 
v, so that E = 0. This result is not surprising if E is interpreted as a measure of non-linear 
gravitational energy, since plane waves propagate linearly. 

Colliding plane wave spacetimes, introduced by Szekeres (1970, 1972) and reviewed 
by Griffiths (1991), have non-zero E, with no definite sign. As an example, the solution 
of Khan and Penrose (1971), describing the collision of two impulsive gravitational waves, 
has E < after the collision, with E becoming unbounded at the singularity formed by 
the collision. (The exact expression is long and unilluminating.) The usual interpretation 
is that the two incoming waves interact non-linearly, producing scattered radiation ^ 2 
whose magnitude grows without bound, a property dependent on the plane symmetry. The 
energy E provides a measure of the non-linear interaction, being zero before the collision 
and becoming unbounded at the singularity. This illustrates both that E is not a preserved 
quantity, and that it is sensitive to gravitational interactions in vacuum. Similar results are 
obtained for generic colliding plane waves, and for asymptotically plane waves (Hayward 
1990, 1992): E equals up to a factor, with ^2 generically becoming unbounded after 
the collision. 

9. Energy inequalities 

On the question of whether a quasi-local energy should be manifestly non-negative, it 
should be noted that negative-energy matter can be described quite consistently in general 
relativity, and in such a situation it would be preferable for the quasi-local energy also to 
become negative. The examples of §8 show that E is positive for certain physically familiar 
gravitational fields, but that negative E may occur in other circumstances, even in vacuum. 
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From purely quasi-local arguments, the fundamental forms are freely specifiable on any 
one 2-surface, and so it is easy to construct 2-surfaces with negative E. Nevertheless, 
it may be possible to find a positivity theorem based on global assumptions, such as 
asymptotic flatness, global hyperbolicity etc, together with a local energy condition on the 
matter. A weaker possibility would be positivity for spacetimes sufficiently close to flat 
spacetime (Christodoulou and Klainerman 1989), or for sufficiently round spheres {e.g. 
Christodoulou and Yau 1988). Conversely, E may be a suitable quantity to control the 
gravitational field in the context of global existence theorems or cosmic censorship. In this 
context it is interesting to note that in the 2-dimensional dilaton gravity theory of Callan 
et al (1992), positive gravitational energy is associated with trapped spatial singularities, 
and negative energy with naked singularities (Hayward 1993b). In spherical symmetry, 
there is a similar relationship between the sign of E and the signature and trapping of 
singularities (cf Christodoulou 1991). 

A related question is that of monotonicity, e.g. whether E increases as A increases 
along a foliation of 2-surfaces. Here it may be useful to consider well- oriented surfaces 
such that 99 < 0, which have an orientation determined according to which expansion is 
negative and which positive. For such surfaces, the Hawking mass decreases internally and 
increases externally (Eardley 1979), which can be interpreted as a quasi-local generalisation 
of the Bondi mass-loss result. A similar result for E would be of interest, but seems more 
difficult to obtain. 

Similar comments apply to energy inequalities such as the isoperimetric inequality 
(Gibbons 1972, Penrose 1973) and the hoop conjecture (Thorne 1972) for trapped surfaces. 
For a marginally trapped sphere, 99 = and (7 = 0, 



and consequently the isoperimetric inequality 16nE 2 > A can be violated by a suitable 
quasi-local choice of shears and twist. Again, global assumptions may be relevant. Note 
however that the twist and internal shear vanish on a Killing horizon (Tod 1986), so that 
a static black hole satisfies the isoperimetric equality 16nE 2 = A on the horizon, a more 
precise result than hitherto obtained (Ludvigsen and Vickers 1983b, Tod 1985, 1986, 1992). 
The very triviality of this result suggests that E may be an appropriate definition for such 
purposes. Expressed equivalently, E is the irreducible mass on a static horizon, i.e. the 
mass of a Schwarzschild black hole with the same area (Christodoulou and Ruffini 1971). 



As noted in the Introduction, there is now a plethora of suggested quasi-local energies, 
which disagree even in the simple Reissner-Nordstrom and Kerr cases (Bergqvist 1992a). 
The ambiguity can only be resolved by settling on generally agreed criteria, and I would 
like to make a few such suggestions. Most fundamentally, the energy should be well-defined 
on any compact orientable spatial 2-surface, or at least for the zero genus case. This firm 
but fair requirement disqualifies almost all of the contenders, which depend on particular 




10. Criteria for quasi-local energy 
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symmetries (Komar 1959, Katz et al 1988, Kulkarni et al 1988, Bergqvist and Ludvigsen 
1989), a preferred 3-surface (Geroch 1973, Nester 1991, Brown and York 1993), global 
requirements (Ludvigsen and Vickers 1983a, Bartnik 1989) or the existence of solutions to 
particular equations (Penrose 1982, Dougan and Mason 1991, Bergqvist 1992b). To the 
best of my knowledge, this leaves only the Hawking mass as an unambiguous quasi-local 
definition. 

The other basic criteria are that the energy vanishes in flat spacetime, and yields the 
Bondi mass at null infinity and the adm mass at spatial infinity. It would also be of interest 
to investigate the Newtonian and linearised limits. Eardley (1979) and Christodoulou and 
Yau (1988) give a comparable list of criteria, of which positivity and monotonicity remain 
as open questions for E. Bergqvist (1992b) also gives a list of criteria, and shows that 
there are infinitely many possibilities which satisfy it. Such extravagance can be curtailed 
by demanding that the definition should depend on the fundamental forms only, as these 
are the only purely geometrical quantities associated with a 2-surface. The fundamental 
forms are also the dynamically independent parts of the gravitational field on a 2-surface, 
or equivalently the free gravitational data. 

Even with such criteria, the arguments of §7 show that a unique 'correct' quasi-local 
energy will require further (presumably dynamical) justification. Some of the suggested 
quasi-local energies involve spurious dynamics, in the sense that they are based on the 
introduction of fields satisfying various particular equations other than the Einstein equa- 
tions. The main point of this article is that the '2+2' Hamiltonian 7i is a dynamically 
preferred quantity which can justifiably be interpreted as the energy density of the Einstein 
gravitational field referred to a 2-surface element. The fact that the associated quasi-local 
energy E has the agreeable properties described herein may be taken as firm supporting 
evidence. In view of this, are there any objections to using the '2+2' Hamiltonian of a 
field to define its quasi-local energy? Is E the long-sought gravitational energy? 
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Appendix: the Bondi mass 

When defining the Bondi mass as a limit of a quasi-local energy along a null foliation 
of 2-surfaces S, there is an important subtlety, namely that the limit depends on the 
foliation used, and in particular whether the foliation is based on an affine parameter or 
a luminosity parameter (area radius). Conversely, different quasi-local integrals must be 
used in order to obtain the same limit. f The following is based on a calculation of Jorg 
Frauendiener (private communication). 

Assume a parameter s which is related to the area radius r by 

r = s + y S - 1 + 0(s- 2 ) 

for some /. Take a spin-basis (l,n,m,fh) such that m spans S and / = d/ds. The 
Newman-Penrose convergence p is then given by p = — s _1 + fs~ 3 + 0(s -4 ). 

The original definition of the Bondi mass (Bondi et al 1962, Sachs 1962) involved 
the Bondi coordinates (w, r, <p), where u is null and the area 2-form of a constant-(-u, r) 
surface is p = r 2 sin-tfcM A dp. Namely, the definition is 

M B = — I MsintfcM A d<p, 



4n _ 

where the mass aspect M is defined by the expansion 

— = -1 + 2Mr _1 + 0(r" 2 ), 



or equivalent ly by 



M = i lim r hi 

' gVU 



2 

i — >oo 



-.rr 



This may be translated in terms of the parameter s by 

g rr = g ss {l _ l /s - 2) 2 + g su fuS -l + g ru = g su {1 _ l /fl - 2) + 

Consider two special choices of s. Firstly, take s to be the affine parameter f used in §9.8 
of Penrose and Rindler (1986), for which 

/ = -a ^°, /" = 1, = -1 -2m° 2 r~ 1 +0(f" 2 ), 

where 9ft denotes the real part, 

a = a f- 2 + O(f- 3 ), ^ 2 = ^f- 3 + 0(f" 4 ), 



f I am grateful to Jorg Frauendiener and John Stewart for clarifying an earlier confusion 
on this point. 
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and a and ^2 are the Newman-Penrose shear and complex 'Coulomb' term respectively. 
Hence the mass aspect is given by 

M = -m° 2 + \f u = + &°^ u ) = -M(v° 2 - &°& '). 

Alternatively, taking s to be the area radius r, 

/ = 0, 9 — = -l-2mlr- 1 +0{r- 2 ), 

and so simply 

m = -m%. 

There are various equivalent expressions (e.g. Ashtekar and Magnon-Ashtekar 1979), since 
h ac h bd C abcd = 2e- m h ac l b n d C abcd = -2 e - 2m l a n b l c n d C abcd = -4m 2 , 

using the null normals l a and n a (§2). 

The above calculation shows that the 'Coulomb' terms for the spin-bases adapted to 
r and f are related by 

m° 2 = K(^-a°a 0/ ), 

essentially because 

r = r-^a ^°f- 1 +O(f- 2 ). 
Consequently, the Bondi mass can be expressed in either of the forms 

M B = - lim TTVT^r / ^*2 = -lim 7TV TF - / A^(*2 - ^a')- 

r^oo 2w V lD7T J s r^oo 2w V 1D7T J§ 

These are limits of the quasi-local energy E and the Hawking mass Mh respectively, 

M B = lim E(S) = lim M H (S). 

Thus the Bondi mass is the limit of the Hawking mass along the affine foliation used 
by Penrose and Rindler (1986), and is also the limit of the quasi-local energy E along a 
foliation parametrised by area radius. The area radius or luminosity parameter was the 
original choice (Bondi et al 1962) and an affine parameter is used in most other work 
(e.g. Newman and Unti 1962, Newman and Penrose 1968, Dixon 1970, Newman and Tod 
1980, Penrose and Rindler 1986, Stewart 1990). The area radius is in many ways the most 
natural choice, for instance from purely geometrical considerations, or from the link with 
the adm mass (§6). Additionally, note that other affine parameters would also require 
different expressions, with the particular affine parameter f being chosen to agree with the 
area radius to the highest possible order. 
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